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Abstract 

^-nucleoli interactions are investigated in the framework of coupled-channels 
dynamics based on the next-to-leading order chiral SU(3) meson-baryon ef- 
fective Lagrangian. A recent determination of the Is shift and width of kaonic 
hydrogen enables us to set accurate constraints on the coupled-channels 
meson-baryon amplitudes in the strangeness 5* = — 1 sector. Theoretical 
uncertainties in the subthreshold extrapolation of the coupled-channels am- 
plitudes are discussed. Using this framework, we give predictions for K~- 
neutron interactions and for the spectrum of the A(1405) resonance. A 
simplified, effective three-channel model using leading order chiral SU(3) 
meson-baryon interactions is also constructed for convenient application in 
-S'-nuclear few-body calculations. 
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1. Introduction 

An important and challenging theme in strangeness nuclear physics is the 
dynamics of antikaons interacting with nucleons and nuclei. The KN interac- 
tion at low energy is strongly attractive and generates the A(1405) resonance 
as a quasi-bound state embedded in the vrS continuum below KN thresh- 
old [T] . One therefore expects that interesting phenomena will also take place 
when the antikaon is injected or stopped in nuclear systems. Much work has 
recently been devoted to investigations of K few-nucleon systems |2l |3l H] and 
of possible bound states of a ^ in heavier nuclei [5]. So far, however, the un- 
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certainties in the subthreshold extrapolation of KN interactions, apart from 
experimental ambiguities, have prohibited firm and consistent conclusions. 

The basis for these studies is the KN two-body interaction. In addition to 
the strong attraction in the elastic KN channel, the other prominent feature 
is the almost equally strong coupling in the transition amplitude for -H- 
ttS. Such strong coupled-channels dynamics is well treated by a unitary 
approach starting from chiral SU(3)ii;xSU(3)L effective field theory [HI [71 El 

Elllo]. 

In earlier work, the data base for low-energy KN interactions came from 
scattering experiments performed from the 1960s to the 1980s. However, the 
K^p total cross section data in elastic and inelastic channels and threshold 
branching ratios did not constrain the scattering amplitude sufficiently well, 
so that the extrapolation of the KN amplitude to the subthreshold energy 
region still suffered from large uncertainties [8^ |9j. 

Information about KN threshold physics comes primarily from measure- 
ments of the energy shift, and width, F, of the Is state in kaonic hy- 
drogen. From these measurements, the real and imaginary parts of the K~p 
scattering length, a{K~p), can be deduced applying the improved Deser- 
Trueman formula [11] (see also Ref. for a precise calculation of the atomic 
level shifts). A precise determination of a{K~p) is crucial for a reliable sub- 
threshold extrapolation of the corresponding scattering amplitudes. Several 
previous kaonic hydrogen measurements [131 E] extracted AE and F, but 
with still large uncertainties. Moreover, a possible quantitative inconsis- 
tency between the DEAR measurements [H] and the scattering data was 
pointed out in several theoretical papers [H HI |TTj. The situation has now 
improved significantly with the advent of the new SIDDHARTA measure- 
ments of kaonic hydrogen [12]. Not only are these data far more precise than 
the previous ones, but as reported in very recent theoretical work [16], they 
are also found to be fully consistent with the existing scattering data. This 
new result enables us to update the theoretical description of the KN in- 
teraction and to reduce the uncertainties in the subthreshold extrapolation 
of the KN amplitude. At the same time it is possible to predict improved 
values for the K~ neutron scattering length, a{K~n), and to set constraints 
for the i^'~-deuteron scattering length [17i| . 

In this paper, we present a further extended study of the KN-ttTj in- 
teraction using the chiral coupled- channels framework with constraints from 
the SIDDHARTA measurement together with the existing scattering data. 
The interaction kernel is constructed based on input from the chiral SU(3) 
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meson-baryon effective Lagrangian. In addition to the results already re- 
ported in our previous letter [16], we discuss implications for the structure 
of the A(1405) and K~n scattering. 



2. Chiral SU(3) dynamics at next-to-leading order 

2.1. Chiral effective Lagrangian 

Our starting point is the chiral SU(3)ijxSU(3)2, meson-baryon effective 
Lagrangian at next-to- leading order (NLO): 

C,ff{BM) = Cm{U) + 4/k^,W) + , (1) 

where Cm^M) with hi = = exp[iA/2$//] is the non-linear chiral meson 
Lagrangian incorporating the octet of pseudoscalar Nambu-Goldstone bosons 
{tt,K,K,t]) in the standard 3x3 matrix representation $. At this stage 
/ ^ 86 MeV is the pseudoscalar decay constant in the chiral limit. 

The baryon octet fields (A^, A, E, S) are collected in the 3x3 matrix B. 
The pseudoscalar meson octet couples to the baryons through the mesonic 
vector current 

v^' = ^^(u^d^u + ud^u^) , (2) 

and the axial vector current 

= j^iu^d^u - = -^u^' . (3) 

The most general form of the meson-baryon interaction Lagrangian at leading 
order 0{p) in the chiral expansion is given by 

4i)^ = Tr(i3(i7^P'^-Mo)S-Di37^75{a^i3}-Fg7M75K,S]) , (4) 

with the chiral covariant derivative V^B = d^B + i[f'^,;B]. Here D and F 
are the low energy constants of the axial vector couplings and Mq is the 
baryon octet mass in the chiral limit. The vector current coupling to the 
baryons involves even numbers of pseudoscalar mesons, while the axial vector 
vertices involve odd numbers of mesons. At next-to-leading order, 0{p'^), the 
Lagrangian introduces several low-energy constants {bi and dj) as 

C^Sb =boTT{BB)TT{x+)+bnTT{B{x+,B})+bpTT{B[x+,B]) 
+ d^ Tt{B K, K, 13]}) + d2 Tt{B [u^, K, B]]) 
+ 4 Tr {B Uff) Tr {B u") + d^ Tr {B B) Tr {u^ u^") , (5) 



3 




Figure 1: Feynman diagrams for the meson-baryon interaction: Tomozawa- Weinberg term 
(i), direct and crossed Born terms (ii) and (iii), and NLO terms (iv). Dashed (soUd) hnes 
represent the pseudoscalar octet mesons (octet baryons). 



where 

,, = _!M(„M. + .tM.t) (6) 

is the exphcit symmetry breaking term with the chiral condensate (0|gg|0) 
and the quark mass matrix JH = diag(m„, m^, m^). 

At tree level in chiral perturbation theory (ChPT), the low-energy con- 
stants bo, })£, and hp are related to the baryon octet masses. The present 
analysis goes systematically beyond tree level and utilizes chiral SU(3) effec- 
tive field theory to determine the interaction kernel of the coupled-channels 
scattering equations in which this matrix kernel is iterated to all orders. 
Therefore the low energy constants need not be identical to those in ChPT 
once the renormalization of the one-particle irreducible graphs is properly 
taken into account. 



2.2. Meson-baryon interactions 

Consider meson-baryon scattering using the chiral Lagrangians Q and 
(jsj). In the SU(3) sector, several meson-baryon channels are coupled in sectors 
with given quantum numbers. The interaction matrix elements are written as 
Vij = {i\V\j) with final and initial channel indices, i and j. The strangeness 
S = —1 and charge Q = sector involves ten channels labeled by the indices 
i = 1,...,10 in the order K'p, K°n, 7r°A, 7r°S°, vr+S", 7r-S+, r/A, ryE", 
K^'E~, K^S^. The interaction Vij is a function of the meson-baryon center- 
of-mass energy a/s, the scattering angles Q = {0,ip}, and the baryon spin 
degrees of freedom cTj and aj. Projecting onto s-wave, the interactions depend 
only on a/s: 

Vij{^) = I dnVij{^s,n,ai,aj) , (7) 
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where the spin average over a-i = aj = ±1/2 is taken. 

The covariant derivative in Eq. Q generates the Tomozawa- Weinberg 
contact term shown in Fig. [l](i). This is the leading-order contribution to 
pseudoscalar meson-baryon scattering, given a^ 

^(TW) 

V^™(v^) = -^MA/-,(2v/^-M,-M,) , (8) 

where /j, Mj and Ei = ^/Mf~+~qf are the meson decay constant, the (phys- 
ical) baryon mass and the baryon energy in channel i. The center-of-mass 
three-momentum qi in that channel is: 



2^s ' ^> 

with the meson mass mj. The normalization factor is Mi = y/Mi + Ei. The 
constants C^J"^^ are determined by SU(3) Clebsch-Gordan coefficients and 
given in Refs. O |8] . 

Next in the ChPT hierarchy are the Born terms derived from the meson- 
baryon Yukawa vertices in Eq. Q. The direct Born term shown in Fig. [l|ii) 
is given by 



k=l 



with channel k denoting an intermediate baryon. The coefficients C*-^j^™°^ 
include the constants D and F and are tabulated in Ref. [8]. The crossed 
Born term (Fig. [l|^iii)) involves the u-channel Mandelstam variable u and is 



^The normalization convention used here is the same as in Ref. [8], with dimensionless 
Vij . It differs from the one used in Ref. [3 [10] by a factor ^jMiMj. 
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given by 



^{Bom)^{Born) 



k=l 
X 



12/./, 



u 



Ml 



M, 



{Mi + Mk){M, + Mk 



s - Mfc + Mi + M. 



s + Ml 



2EiEj 



X In 



(y/^ - Mk + Mi + Mj)} 

s + Ml- mf - m] - 2EiEj - 2^^^/ 



2NfNf 



s + Ml 



2 



2EiEj + 2qiqj 



(11) 



Finally we turn to the next-to-leading oder terms (Fig. [T|iv)). The cor- 
responding interactions are derived from the four-point vertices in Eq. ([5]): 



T/(NLO) 
ij 



fifj 



^(NLOl) ^^(NL02) 



EiEi + 



(12) 



The coefficients (^j]^^*-*^) and Q^^^'^'^^ 



are again summarized in Ref . [8j . They 
include the NLO low-energy constants b^, b^, bp and di. Since we take into 
account the renormalized (physical) masses of the baryons, these low-energy 
constants are not identical to the ones from tree-level ChPT. 



2.3. Chiral coupled- channels dynamics 

In contrast to the SU(2) pion-nucleon systems close to threshold, the 
meson-baryon channels in the strangeness sector are strongly interacting. In 
particular, the KN interaction is sufficiently strong to produce the A (1405) 
as a quasi-bound state below the KN threshold. In such a situation, a 
perturbative calculation does not work and a nonperturbative resummation 
is mandatory in order to account for the strong coupled-channels dynamics. 
The meson-baryon interactions derived above are thus used as the interaction 
matrix kernel V of the coupled-channels Bethe-Salpeter equation for the T- 
matrix: 

T = V + V-G-T = (V"^-G)"^ . (13) 
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Here G is a diagonal matrix, Gij = Gi{Q) Sij, with elements Gi representing 
the meson-baryon loop function in channel i: 



(27r)4 [(g - ky - M2 + ie] (P - mj + ie) ' 
Its logarithmic divergence can be tamed by dimensional regularization as 



(14) 



327r2 



167r2 



s 



1 + ^ artanh , , , ^ , ^ 
yTs \im, + M,)2 - s 



(15) 



The subtraction constants ai(/i) act as renormalization parameters at a scale 



fx such that the T-matrix (13) is scale independent. Note that the last 



expression in Eq. (13) is well justified in formal scattering theory, e.g. by 



using the N/D method which automatically guarantees the unitarity of the 
scattering amplitude. As demonstrated in Refs. [TJ the interaction kernel 
V can be identified with the one derived from tree-level ChPT up to next- 
to-leading order with no problem of double counting. 

2.4- Observables 

The meson-baryon coupled- channels T-matrix is the starting point for the 
calculation of various K~p scattering observables. This T-matrix is related 
to the forward scattering amplitudes fij as 



fijiVs) 



in 



(16) 



The K~p elastic scattering amplitude at threshold defines the scattering 
length, a{K~p) = /ii(a/s = nix- + Mp), a complex number because of the 
absorptive channels converting K~p into ttS and ttA. The energy shift and 
width of the Is state of kaonic hydrogen are related to the K~p scattering 
length, with important second order corrections, as follows fTT]: 



(17) 



AE - ir/2 = -2a^ /i^ a{K'p) [l + 2a /i^ (1 - In a) a{K'p 



where a is the fine-structure constant and the K p reduced mass is given by 
fir = rriK-Mp/ [rriK- + Mp). 
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The total reaction cross sections in the various meson-baryon scattering 
channels are given by 



where the cross section is defined for y/s > Mi + rrii, above the threshold 
of the final-state channel i. For the K~p elastic cross section, we also take 
into account electromagnetic interactions which are important near the K~p 
threshold [S]. The Coulomb interaction gives an additional contribution to 
the diagonal amplitude in the K~p channel: 

fCoul ^ 



2gf as sin2(^em/2) 

T{l-i/{q,aB)) / 2i , .9, 



r(l + i/(giaij)) \qiaB 



expf In sin , (19) 

Voi Or 2 / 



with as = 84 fm, the Bohr radius of the K~p system, and 6cm denoting 
scattering angle. This Coulomb amplitude is added to the strong interaction 
amplitude and the scattering angle is integrated up to cos6'cm < 0.966 to 
avoid the divergence at 6'cm = 0. 

Several combinations of K~p inelastic yields at threshold are known in 
the form of branching ratios defined as 

_ T{K-p^n+^-) _ ^51 p _ T{K-p^7T^A) a^i 



T{K p — TT 11+) (Tei ^{K p — j- neutral states) a-n + a^i 

T{K~p — )■ all inelastic channels) a^i + cr4i + (Xsi + a^i ' 
with all partial cross sections cTij calculated at the K~p threshold. 



[eV] 


r[eV] 


7 


Rn 


Rc 


283 ± 36 ± 6 


541 ± 89 ± 22 


2.36 ±0.04 


0.189 ±0.015 


0.664 ±0.011 



Table 1: Experimental observations of the energy shift and width of the Is state of kaonic 
hydrogen {AE and F) il5i, threshold branching ratios (7, i?„ and Rc) [TS] . 
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3. Results and discussion 



3. 1 . Fitting procedure 

We now describe the systematic fitting procedure used in tlie framework 
of the chiral SU(3) dynamics at NLO level. We first summarize the empir- 
ical constraints that enter this study. Important constraints are the kaonic 
hydrogen shift and width from the SIDDHARTA measurements [15J: 

= 283 ± 36(stat) ± Q{syst) eV , T = 541 ± 89(stat) ± 22{syst) eV . 



The threshold branching ratios ( 20 ) are determined [TH] from K~ capture on 
hydrogen as listed in Table [T] In addition to these K~p threshold constraints, 
we also make use of the total cross section data accumulated in Refs. [19] 
(see Fig. [2]). In contrast to several previous studies, the currently available 
vrE mass spectra are not included in the empirical data base. At this point 
a meaningful comparison with experimental data would require a detailed 
investigation of the different reaction mechanisms generating such spectra 
in each given experiment [in]. Instead, the vrE spectrum emerges here as a 
prediction following the fitting to the previously mentioned quantities. 

The fits to the data base have been performed using three consecu- 
tive schemes, systematically improving the interaction kernel in each step. 
The first setup involves just the Tomozawa- Weinberg (TW) term (|8]), the 
dominant component of the leading-order interaction. In the second step we 



include in additon the direct and crossed Born terms, Eqs. (10) and (11), 
completing the 0{p) ChPT interaction (TWB). The third step incorporates 
all terms (12) of the full NLO model. We define the overall following 



Ref. [9]. Consider a measurement, labeled and the number of data points 
Ui associated with this measurement. We first calculate Xi for each i-th mea- 
surement, then collect all xj multiplied by proper weight factors, and finally 
obtain overall x^/d.o.f., 

xVd.o.f.=— ^^x:^' (21) 

where R is the total number of measurements considered. Furthermore, p 
is the number of parameters that appears in each step, increasing from TW 
via TWB to the full NLO setup. If all measurements have same numbers of 



data points, Eq. (21) reduces to the standard definition of x /d.o.f. 
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In the numerical calculations we use the physical masses for mesons and 
baryons. This is necessary in order to reproduce the correct threshold ener- 
gies in the different meson-baryon channels. The constants associated with 
the axial vector baryon couplings are given hj D = 0.80 and F = 0.46 (i.e. 

= D + F = 1.26). The meson decay constants are chosen at their physical 
values 



U = 92.4 MeV , /x = (1.19 ± 0.01) , f^ = (1.30 ± 0.05) U . (22) 

Using the physical masses and decay constants is justified by taking into 
account the renormalization of the baryon masses and meson fields. As a 
consequence, parts of the effects of the NLO parameters Bq, ho and hp have 
already been absorbed in the renormalized quantities, shifting the baryon 
octet masses from their degenerate chiral limit, Mq, to their physical values; 
the parameters used in the interaction kernel are therefore to be interpreted 
as the renormalized ones. We denote these renormalized NLO parameters as 
ho, ho and hp- They are expected to be considerably smaller in magnitude 
than the ones usually quoted in tree-level chiral perturbation theory. 
Thus, the free parameters to be determined by the fits are: 

i) the subtraction constants aj(/i) that are assumed to be isospin symmetric; 

ii) the low energy constants in the NLO Lagrangian h^, bp,, hp and di. 
Note that the set i) appears in all three TW, TWB and NLO versions, 
whereas set ii) is only used in the full NLO scheme. In the fitting proce- 
dure we also allow for small modifications of the K and r] decay constants 



within the uncertainties given in Eq. (22). 



3.2. Fit results 

With the TW terms alone a reasonable overall fit can already be reached 
with x^/d.o.f. = 1.12, as shown in Table [2] However, although the branch- 
ing ratios and cross sections are reproduced quite well as shown in Table [3] 
and Fig. |2| the kaonic hydrogen energy shift {/S.E = 373 eV) exceeds the 
empirical bound. It should also be noted that the fitted values of the sub- 
traction constants become large in the ttA and r^E channels, beyond their 
expected "natural" size, \a\ ~ 10~^ [3 El]. This result indicates the limit of 
applicability of the simple TW model for the quantitative discussion of the 
scattering amplitudes, given the remarkable accuracy that has been achieved 
in the SIDDHARTA measurements. 

Next we consider the TWB scheme in which the direct and crossed Born 
terms are combined with the TW term. In this case, the x^/d.o.f. changes 
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TW 


TWB 


NLO 


axN (10-3) 




-1.57 


-1.04 


-2.38 


a^A (10-=^) 




-107.97 


-8.06 


-16.57 


a.E (10-3) 




2.31 


2.96 


4.35 


ar,A (10-3) 




-0.20 


-3.46 


-0.01 


a^E (10-3) 




216.37 


3.52 


1.90 


axs (10-3) 




39.48 


12.51 


15.83 


/x (MeV) 




110.8 


109.0 


110.0 


/, (MeV) 




124.5 


124.6 


118.8 


bo (10-2 GcV- 


^) 






-4.79 


fez) (10-2 GcV- 


-^) 






0.48 


bp (10-2 GeV- 








4.01 


di (10-2 GeV- 








8.65 


c^2 (10-2 GcV- 








-10.62 


^3 (10-2 GcV- 








9.22 


(^4 (10-2 GeV- 








6.40 


xVd-O.f. 


1.12 


1.15 


0.96 



Table 2: Parameters resulting from the systematic analysis, using leading order (TW) 
plus Born terms (TWB) and full NLO schemes. Shown arc the isospin symmetric subtrac- 
tion constants aj(/x) at /z = 1 GeV, the meson decay constants Jk and /,,, the renormalized 
NLO constants bi and dj, and x^/d.o.f. of the fit. 
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Figure 2: Calculated K~p elastic, charge exchange and strangeness exchange cross sec- 
tions as functions of laboratory momentum, compared with experimental data |19| . 
The dashed, dotted and solid curves represent TW, TWB and best fits of the full NLO 
calculations, respectively. 
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TW TWB NLO 



AE [eV] 
r [eV] 



373 
495 



377 
514 



306 
591 
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2.36 
0.20 
0.66 



2.36 
0.19 
0.66 



2.37 
0.19 
0.66 



pole positions 
[MeV] 



1422 - 16 i 
1384 - 90 i 



1421 - 17i 
1385 - 105 i 



1424 - 26 i 
1381 -81i 



Table 3: Results of the systematic analysis using leading order (TW) plus Born terms 
(TWB) and full NLO schemes. Shown are the energy shift and width of the Is state of 
kaonic hydrogen (AE and F), threshold branching ratios (7, i?„ and Rc), and the pole 
positions of the isospin / = amplitude in the KN-ttY. domain. 



only marginally. The energy shift AE = 377 eV is not improved by the 
inclusion of the Born terms, but the values of the subtraction constants 
(Table |2]) are now reaching "natural" sizes. From the theoretical point of 
view, this now indicates a consistent description of the interaction kernel 
and the loop function. 

Finally the NLO terms are added in the construction of the full ampli- 
tudes. The analysis provides the best fit parameters with x^/d.o.f. = 0.96. 
The kaonic hydrogen shift and width are now within the error bars of the 
SIDDHARTA measurements. Threshold branching ratios and total cross 
sections of the K~p scattering and reaction processes are well reproduced as 
demonstrated in Table [3] and Fig. |2} The parameters determined by the NLO 
fit are altogether meaningful: natural-sized subtraction constants and small 
renormalized NLO parameters. The stepwise improvement of the theoretical 
description of the meson-baryon amplitudes in the three schemes from TW 
via TWB to NLO is evident, emphasizing the important role of the accurate 
kaonic hydrogen data in constraining chiral SU(3) dynamics. In contrast, 
the scattering data alone (see Fig. [2]) do not provide a sensitive test for the 
different schemes. 
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To estimate the uncertainty in the subthreshold extrapolation, we ex- 
amine variations of the parameters around their best fit values. We vary 
the subtraction constants with the condition that the shift and width of the 
kaonic hydrogen are reproduced within experimental errors. In addition, the 
total cross section of the K~p — vr^A process is also used to set constraints 
on the 1 = 1 amplitudes. With these error assignments, the allowed ranges 
for the subtraction constants are taken at the scale /x = 1 GeV are found 
to be: 

n- — 00+0.17 ^ , _ lfiP;7+16.81 ^ _/| qt;+0.31 

a^A = - 0.0i;?:22 , a,s=1.90l°J^ , a^^ =15.831};^^ , (23) 

in units of 10^'^. The corresponding error bands in the cross sections are 
shown in Fig. |3} together with the best-fit curves. It is evident that the cross 
sections are well constrained by the kaonic hydrogen data: there is mutual 
consistency between scattering and threshold measurements. 

3.3. Predictions from improved chiral SU(3) dynamics 
3.3.1. K^p subthreshold amplitudes and structure o/A(1405) 

Once this consistent theoretical description of all K^p observables is ob- 
tained, it is instructive to perform extrapolations of the amplitudes to sub- 
threshold and complex energies. This is particularly important for the un- 
derstanding of the A(1405) resonance as a quasibound K~p (/ = 0) state 
embedded in the ttE continuum, as well as for the far-subthreshold KN in- 
teraction that is relevant in the context of possible ^-nuclear clusters. In 
Fig. |4], we show the subthreshold extrapolation of the real and imaginary 
parts of the K~p K~p amplitude by the best-fit NLO scheme and the 



uncertainty bands constrained by Eq. (23). As seen in the figure, the am- 
plitude exhibits the structure of the A(1405) resonance emerging from the 
strong attracttion in the J = component of the amplitude. We find that 
the subthreshold extrapolation is stable thanks to the accurate constraint at 
threshold. 

3.3.2. The two-poles scenario 

Next we look for poles in the second Riemann sheet of the complex energy 
plane to study the coupled-channels structure of the A(1405) resonance. With 
the best-fit result in the NLO scheme, pole singularities between the KN and 
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Figure 3: Calculated K~p elastic, charge exchange and strangeness exchange cross sections 
as functions of K~ laboratory momentum, compared with experimental data |19| . The 
solid curves represent best fits of the full NLO calculations to the complete data base 
including threshold observables. The shaded uncertainty bands are explained in the text. 
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7rS thresholds are found at 



^1 =1424-126 MeV , 22 =1381 - i 81 MeV . 

The higher energy Zi pole is dominated by the KN channel and the lower 
energy Z2 pole receives stronger weight from the vrS channel. This confirms 
the two-poles scenario of the A(1405) [71 |22l [23]. Actually, the existence of 
two poles around the A (1405) resonance had been found in previous NLO 
calculations [H [9], but the precise location of the poles, especially of the 
lower one, could not be determined in these earlier studies, given the lack of 
precision in the empirical constraints. 

In the present analysis, the SIDDHARTA measurement provides much 
more severe constraints also on the pole positions. The real parts of zi and 
Z2 are remarkably stable in all three TW, TWB and NLO schemes. The 
imaginary parts deviate within < 20 MeV between these schemes, as seen 



in Table |3} Using the error analysis from Eq. (23) together with the best-fit 
NLO results, one finds: 

^1 =1424t^3 -i26t?4 MeV , ^2 =13811^^ - i 8ll^^ MeV . (24) 

The uncertainties of the pole locations are thus significantly reduced from 
previous work, and the two-poles structure of the A(1405) is now consistently 
established with the constraints from the precise kaonic hydrogen measure- 
ment. Because of isospin symmetry, the two poles are stable against varia- 
tions of the 1=1 subtraction constants (the ones in the ttA and r/S channels). 
The error assignments in the pole positions and half widths are mainly re- 
flecting the uncertainties of the KN and vrS subtraction constants. 

3.3.3. K~p and K~n scattering lengths 

A discussion of low-energy ^-nuclear interactions requires the knowledge 
of both the K~p and K~n amplitudes near threshold. The complete KN 
threshold information involves both isospin / = and 1 = 1 channels. The 
K~p scattering length a{K~p) = [a^ + ai]/2 is given by the average of the 
/ = and 1 = 1 components, whereas the K~n scattering length a{K~n) = 
tti is purely in / = 1. Note that Coulomb corrections to a{K~p) and isospin 
breaking effects in threshold energies may be significant [H] and must be 
taken into account in a detailed quantitative analysis. 

We first extract the scattering length a{K~p) from the SIDDHARTA 



measurements [15] using Eq. (17). The result is: 



Re a(i\:-j9) = -0.65 ±0.10 fm , lma{K-p) = 0.81 ±0.15 fm , (25) 
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1340 1360 1380 1400 1420 1440 

Vi [MeV] 



1340 1360 1380 1400 1420 1440 
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Figure 4: Real part (left) and imaginary part (right) of the K~p — > K~p forward scatter- 
ing amplitude obtained from the NLO calculation and extrapolated to the subthreshold 
region. The empirical real and imaginary parts of the K~p scattering length deduced from 
the recent kaonic hydrogen measurement (SIDDHARTA [T3]) are indicated by the dots 
including statistical and systematic errors. The shaded uncertainty bands are explained 
in the text. 



where the uncertainties reflect the experimental errors. The predictions from 
chiral SU(3) dynamics, proceeding again through the sequence of TW, TWB 
and full NLO schemes, gives the following values for the K~p scattering 
length: 

a{K-p) = -0.93 + 10.82 fm (TW) , (26) 
a{K-p) = -0.94 + 10.85 fm (TWB) , (27) 
a{K-p) = -0.70 + 10.89 fm (NLO) . (28) 

The large magnitude of Re a{K~p) in the TW and TWB schemes corresponds 
to the overestimation of the kaonic hydrogen energy shift in these approaches, 



while the best-fit NLO result is fully compatible with the value (25) deduced 
from the experimental data. 

To calculate the K~n scattering length, we construct the coupled-channels 
amplitudes in the charge Q = —1 sector {K~n, tt~A, tt^S*^, tt^S", r^S" and 
K^S'), again using physical meson and baryon masses in order to take into 
account isospin breaking effects in the threshold energies. With the same 
subtraction constants as in the Q = sector, the calculated K~n scattering 
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lengths are: 



a{K-n) = 0.29 + i 0.76 fm (TW) , (29) 
a{K-n) = 0.27 + i 0.74 fm (TWB) , (30) 
a{K-n) = 0.57 + i 0.73 fm (NLO) . (31) 

The relatively large jump in Re a{K~n) when passing from "TW" and "TWB" 
to the best-fit "NLO" scheme is strongly correlated to the corresponding 
change in lie a{K~p). Thus, to determine the 1 = 1 component of the KN 
scattering length, it is highly desirable to extract the K~n scattering length, 
e.g. from a precise measurement of kaonic deuterium [2H IT7] . 

Next, consider the subthreshold extrapolation of the complex elastic K~n 
amplitude. Fig. |5] shows the real and imaginary parts of this amplitude. Note 
that the 1 = 1 KN interaction is also attractive but weaker than the / = 
interaction so that f{K~n — )■ K~n) is non-resonant. In the absence of 
empirical threshold constraints for the K~n scattering length one still faces 
relatively large uncertainties. Variation of the subtraction constants within 



the range of Eq. (23) applied to the NLO scheme leads to the following 



estimated uncertainties: 

a(ir-n) = 0.57l°:°f + i0.72;°:26 fm . (32) 

The errors in a{K^n) relate primarily to the uncertainty of the subtraction 
constant in the vrA channel. 




Figure 5: Real part (left) and imaginary part (right) of the K n ^ K n forward 
scattering amplitude extrapolated to the subthreshold region. 
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3.3.4- invariant mass distribution 

One of the important features of chiral SU(3) coupled-channels dynamics 
is the pronounced channel dependence of the A(1405) production spectra 
reflecting the two-poles nature of the A(1405) [22]. To demonstrate this two- 
mode structure, we show the imaginary parts of the scattering amplitudes 
ttS ttS (Fig. g left) and KN KN (Fig. g right) in the / = channel. 
These strength functions exhibit the A(1405) spectrum as seen in different 
channels. Evidently, there is no single universal invariant mass distribution of 
the A(1405). As seen in the figure, the imaginary part of the KN amplitude 
has its maximum close to 1420 MeV, whereas the position of the peak in 
the ttS spectrum is shifted downward from the KN — )■ KN amplitude to 
about 1380 — 1400 MeV. This is a consequence of the strong KN o ttE 
coupled-channels dynamics dictated by chiral SU(3) symmetry. The different 
shapes and positions of the spectral distributions in Fig. |6] represent the 
coupled modes associated with the two poles Zi^2 discussed earlier. While 
the subthreshold KN spectrum has its maximum closer to the location of 
the "upper" pole zi, the ttE spectrum receives a stronger weight from the 
second, "lower" pole Z2. 

The right panel of Fig. |6] includes for reference and orientation the experi- 
mental spectrum of the 7r~S"'" channel in the decay S+(1660) — )■ 7r+(7r~E"'") [23] 
It should however be noted that a direct comparison of this histogram with 
the imaginary part of the calculated / = ttS amplitude is not meaningful. 
The measured spectrum is not pure / = and the relative weights of the 
initial states (vrE, KN , . . .) are not known. In addition, because the energy 
of the three-body 7r+(7r^S+) system is restricted to form the S(1660), the 
higher tail of the vr^S^ spectrum is suppressed because of the small avail- 
able phase space [TD]. It is therefore necessary to construct elaborate reaction 
models to compare the imaginary part of the vrE amplitude with experiments, 
including new vrS spectra recently reported from different experiments |26j . 

3.4. A schematic KN-ttH-tiA model 

For applications in studies of strange dibaryons using sophisticated few- 
body techniques, a tractable model that starts just from the dominant leading 
order Tomozawa- Weinberg interaction terms in a reduced model space of 
KN-ttTi-ttA channels, is often quite useful. The only free parameters in such 
a schematic model are the three subtraction constants, af^^.a-^Y, and a^A, in 
those three channels. Our aim here is to construct an "effective" Tomozawa- 
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Figure 6: Imaginary part of the 7 = KN (left) and ttS (right) amphtudes together 
with error bands permitted by SIDDHARTA experiments. The histogram (arbitrary unit) 
in the right panel denotes the experimental data of the 7r~ E+ spectrum in the decay of 
S+(1660) 7r+(7r-E+) [35]. 

Weinberg (ETW) model that reproduces the results obtained with the best- 
fit, full coupled-channels NLO scheme as well as possible. 

With "natural" -sized isospin-symmetric subtraction constants, a^jy = 
-1.79 X 10~^ a^s = 1-81 x 10"^ and a^A = 7.84 x 10"^ and with the 
meson decay constants = 92.4 MeV and fx = 109.0 MeV, one can indeed 
produce a reasonable set of K threshold quantities as listed in Table |4j The 
double-pole nature of the coupled-channels dynamics is fully maintained in 
this schematic model. The subthreshold K~p — )■ K~p amplitude is shown 
in Fig. [7| in comparison with the best-fit NLO results and their error bands. 
Although some deviations are observed around the KN threshold and the 

cannot compete with the one achieved in the best-fit NLO approach, the 
amplitudes of this ETW model for ^/s < 1425 MeV are well within the 
uncertainties permitted by SIDDHARTA. They compare well with the best- 
fit results in the subthreshold energy region. This simplified ETW model can 
therefore be adopted as input in various practical applications. 

4. Summary 

We have demonstrated that the new kaonic hydrogen measurements, to- 
gether with total cross section data and threshold branching ratios, are suc- 
cessfully described in the framework of chiral SU(3) coupled-channels dy- 
namics with input based on the NLO meson-baryon effective Lagrangian. 
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AE [eV] 


r [eV] 


7 


Rn 


Rc 


pole positions [MeV] 


ETW 


338 


442 


2.26 


0.25 


0.62 


1423 - 22 i 1375 - 65 i 



Table 4: Results from the effective three-channel Tomozawa- Weinberg ("ETW") model. 
Shown are the calculated Is energy shift and width of kaonic hydrogen (A_B and F), 
threshold branching ratios (7, i?„ and Rc), and the pole positions of the isospin / = 
amplitude in the KN-nl] domains. 




1340 1360 1380 1400 1420 1440 1340 1360 1380 1400 1420 1440 



^/5 [MeV] |McV] 

Figure 7: Real part (left) and imaginary part (right) of the K^p — > K^p forward scat- 
tering amplitude together with best-fit results and with error bands permitted by the 
SIDDHARTA measurements. The dashed and solid curves denote the results of "ETW" 
and best-fit NLO, respectively. 

Our systematic study uses physical hadron masses and physical pseudoscalar 
meson decay constants. It reveals the hierarchy of the interaction kernel 
derived in chiral perturbation theory. It is important to point out again 
that the best fit to all existing data has been performed with the constraint 
that the NLO parameters of the chiral SU(3) meson-baryon Lagrangian stay 
within 'natural' limits, such that NLO terms remain small compared to the 
leading-order input. Alternative fits using an unrestricted parameter space 
would be possible, with a x^/d-o.f comparable to our "best fit". However, 
in this case the NLO parameters would turn out unacceptably large and an 
unphysical pole in the 1 = 1 channel below KN threshold would appear as 
a consequence. If this were the only possible option, such a scenario would 
imply that the chiral SU(3) effective field theory coupled-channels approach 
is inconsistent and meaningless for this purpose. The non-trivial observation 
that an optimal fit can be achieved using a consistent hierarchy of LO and 
NLO terms, together with physical values of the pseudoscalar decay con- 
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stants, justifies our conclusion that this can indeed be called a best fit. This 
refined theoretical framework for the KN interaction has several important 
consequences. 

The stringent constraints from the accurate kaonic hydrogen measure- 
ments reduce the uncertainties in the subthreshold extrapolations of the KN 
amplitude significantly. This is an important step towards raising the pre- 
dictive power in calculations of possible dibaryon states in the RNN-ttIIN 
three-body system. The two-poles nature of the A(1405) is confirmed with 
considerably smaller ambiguities than in previous work. The predicted imag- 
inay parts of the KN and ttE amphtudes in the region of the A(1405) show 
a pronounced relative shift in their peak positions, reflecting the two-modes 
scenario of the KN — ttS coupled channels. The K~n scattering length is 
predicted using the same framework. Some ambiguities in the 1 = 1 compo- 
nent of the KN interaction still remain. At present this sector is constrained 
by the poorly known K~p — > 7r°A cross section data and threshold branching 
ratio Rn- The determination of the K'~n scattering length, e.g. through a 
measurement of kaonic deuterium, would be of great importance in order to 
set further constraints in the dynamics of the KN-tt'E system. 
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